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We study the possibility to bind eight quarks in a molecular hadronic system composed of two 
nucleons and an antikaon, with the quantum numbers of a hexaquark flavour, in particular with 
strangeness -1, isospin 1/2, parity -, baryonic number 2 and two possible spins, or 1. These 
exotic hadrons are motivated by the deuteron, a proton-neutron boundstate, and by the model of 
the A(1405) as an antikaon proton boundstate. We discuss the possible production of this hadron 
in the experiments which are presently investigating hot topics like the + pentaquark or the K~ 
deeply bound in nuclei. The K~ • N interactions and the coupling to other channels are computed 
microscopically from a confining and chiral invariant quark model resulting in local plus separable 
Gaussian potentials. The N • N interactions used here are the state of the art Nijmegen potentials. 
The binding energy and the decay rate of the K~ • N and K~ • N • N systems are computed 
with configuration space variational methods. The only systems that bind with our microscopic 
interaction are the K~ • N in the 1 = channel and the K~ • N • N in the S — channel. 



I. INTRODUCTION 

Here we study the binding energy of a possible molec- 
ular antikaon, nucleon and nucleon three-body system. 
Importantly, from a hadronic perspective, a bound K~ • 
N • N can be at most reduced to a hexaquark, a exotic 
system. In a nuclear perspective, the deep binding of an- 
tikaons in nuclei has been predicted and are also actively 
searched. A bound K~ • N • N would constitute a simple 
antikaon- nucleus system, the K _H . 

Moreover, the K~ • N • N can be formed with an- 
tikaon (K~) deuteron (p • n) scattering. Other exotic 
tetraquarks, pentaquarks or hexaquarks are also very 
plausible, but they are all harder to produce experimen- 
tally because they would need at least strangeness and 
charm. The several experiments dedicated to pentaquark 
searches (where not only the Kaon, but also the antikaon 
may interact with nuclei), or to antikaon-nuclear bind- 
ing at RCNP, JLab, KEK, DAFNE and at many other 
laboratories, are already able to search for the proposed 
K~ • N • N '. In particular evidence for K~ • N • N 
has already been found by the FINUDA collaboration 
at DAFNE @, 9. In Fig. Q] different possible produc- 
tion mechanisms are depicted. They are similar to the 
A(1405) production mechanisms, except that the K~ 
scatters on a deuterium nucleus, not on a hydrogen nu- 
cleus. The process in Fig. Q] (a) is only possible if the 
width of the K~ • N • N is oi the order of its binding 
energy. The process in Fig. [1] (b) is always possible, but 
is suppressed by the electromagnetic coupling. Processes 
in Fig. [T](c) , (d) are dominant. The production process 
(c) should occur in experiments designed for the produc- 
tion of the pentaquark + , where the K~ has sufficient 
energy to create a it. The production process (d) should 
occur in atomic Kaon experiments, where the K~ has a 
low energy but where a larger nucleus is used and the 
remaining nucleus may absorb the virtual pion. 

The K~ • N • N binding is suggested by the deuteron 
2 H, a proton-neutron boundstate, and by the model of 
the A(1405) as an antikaon proton boundstate. How- 
ever the subtleties of these two-body subsystems require 
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FIG. 1: Different K • N • N production mechanisms 



a precise computation of the binding energy of the three- 
body system. For instance the binding of the deuteron 
(p • n) requires a d-wave component 4] . Moreover it is 
not clear yet if the K~ • N binds in an / = 1 S baryon. 
Notice that in the 1 = channel there is strong evi- 
dence that the A(1405) is constituted of two poles, re- 
spectively dominated by the K~ • N and the ir • E chan- 
nels @, 0, 1, S M, El, E2, m • 

The K~ • N • N binding is a v ery recent topic. While 
Oset, Toki and collaborators offer a different in- 

terpretation to the FINUDA peak, very recent theoretical 
calculations, fitting the K~»N scattering data with effec- 
tive local potentials [l6|, [13, EBl or with separable poten- 
tials [HI agree with a relatively deep and wide K~»N»N 
resonance. 

Here we specialize in the microscopic quark computa- 
tion of the binding energy and width of the K~ • N • N 
systems. Section II is dedicated to the N»N and K*N in- 
teractions. We compute microscopically the KuN inter- 
action, with a confining and chiral invariant quark model 
and using the Resonating Group Method. In Section 
III we calibrate our interactions with the experimental 
K + • N data, and we study the binding and decay of 
K~ • N and K~ • N • N systems. Finally we conclude. 
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II. FROM QUARKS TO THE K»N 

There are several very precise models of the N • N 
interaction Q , compatible with the experimental N • N 
phase shifts. All contain a long-range attractive poten- 
tial, a medium-range attractive potential, an a hard core 
or short range potential. A picture consistent with QCD 
emerges if the long range attraction is due to the Yukawa 
one-pion-exchange, the medium range is due to two-pion 
or one-sigma exchange, while the short range repulsion 
is dominated by the quark interactions. Different models 
exist because the 3 Si ^ 3 Di coupling is partly equiva- 
lent to the attractive part of the 3 Si potential. Moreover 
the repulsive hard core potential essentially pushes the 
wavefunction outside the repulsive core region, and this 
can be accomplished for different heights of the repulsive 
core. Also, separable or local potentials can be used. 

However the K~ uN interaction has not yet been mod- 
elled with the same detail of the N • N interaction || . 
Notice that not only the Kaons are unstable particles, 
with less experimental results, but also there are large 
inelastic coupled channel effects in the K~ • N scatter- 
ing. Moreover, there is strong evidence that the A(1405) 
does not have a single pole with width T of AOMeV, but 
two narrower poles, one closer to the K~ • N threshold 
and another dominated by the 7r • £ channel. 

This leads us to compute the K~»N interaction micro- 
scopically at the quark level. Here we assume a standard 
Quark Model (QM) Hamiltonian, 



(1) 



where each quark or antiquark has a kinetic energy Tj 
with a constituent quark mass, and the colour depen- 
dent two-body interaction Vij includes the standard QM 
confining term and a hyperfine term, 



= — Xi ■ Xj V con f(r) + Vh yp {r)Si ■ Sj 



(2) 



For the purpose of this paper the details of potential @ 
are unimportant, we only need to estimate its matrix 
elements. The hadron spectrum is compatible with, 
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Moreover we include in the Hamiltonian (TTJ) a quark- 
antiquark annihilation potential Aq. Notice that the 
quark-antiquark annihilation is constrained when the 
quark model produces spontaneous chiral symmetry 
breaking [H E), [H, HI HI [H[ . In the tt Salpeter equa- 
tion, the annihilation potential A cancels most of the 
kinetic energy and confining potential 2T + V, 



(A) s=0 ~(2T+V)s=o-(V hyp } 



(4) 



leading to a massless pion in the chiral limit. We stress 
that the QM of eq. ([T]) not only reproduces the meson 



(EP U ) = Ex 









*> 






(?5 







(Z5 



P). 



{Ms) - 







1 ^ 


s <* 3 




*> 


-1 4 


iA(p-ff\ 






\ ' 



FIG. 2: Examples of our RGM overlaps for the baryon-meson 
effective interaction: (a) norm exchange overlap; (b) kinetic 
exchange overlap; (c) interaction exchange overlap; (d) anni- 
hilation overlap. 



and baryon spectra as quark and antiquark bound-states, 
but it also complies with the PCAC theorems [23l. [241. [2o1. 

[H,[13. 

For the pentaquark system, the Resonating Group 
Method (RGM) [28.J is convenient to arrange the wave 
functions of quarks and antiquarks in antisymmetrized 
overlaps of simple colour singlet quark clusters, the 
baryons and mesons. The effective potential of the 
meson-baryon system is computed with the overlap of 
the inter-cluster microscopic potentials, 

V^a = (<l>B<t>A\-(Vu + V 1 - 5 + 2V24 + 2V 2E )3P 14 

+3A 1 - 5 \<j )A( f )B )/(<t> B (f> A \l - 3P 14 \cf> A q> B ) ,(5) 

where Py stands for the exchange of particle i with par- 
ticle j, see Fig. [5] This results in an algebraic colour x 
spin x flavour factor and a spatial overlap. A convenient 
basis for the spatial wave- functions of the meson A and 
the baryon B is the harmonic oscillator basis, 



0000 (P) =Ma 
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(6) 



In what concerns the quark exchange diagrams, the 
spin independent part of the quark-quark(antiquark) po- 
tential in eq. {2} essentially vanishes because the clus- 
ters are colour singlets. The only potential which may 
contribute is the hyperfine potential. The quark Pauli 
exchange also leads to a positive, repulsive potential. 

Importantly, each annihilation diagram can be related 
to an exchange diagram with the crossing symmetry of 
one of incoming baryon with the outgoing baryon. This 
is evident in Fig. [2](c) and Fig. [2](d). The corresponding 
diagrams are opposite because the annihilation diagram 
does not have the Pauli exchange minus sign of the ex- 
change diagram. Another important difference occurs, 
the exchange diagrams produce non-local simple separa- 
ble potentials, while the annihilation diagrams produce 
simple local Gaussian potentials. Nevertheless all the 
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different diagrams produce hard core potentials propor- 
tional to (Vhyp), repulsive in exchange diagrams and at- 
tractive in annihilation diagrams. Notice that attraction 
may only occur in non-exotic channels. 

We summarize [H 13, iH, [H IH the effective poten- 
tials computed for the relevant channels, 

Vk.n = c K 2 (V hyp )§(l + §T K -r N ) |^oo)(0oool , 
Vk.n(?) = ~c K 2 (V hyp )2V2(l~lr K ■ r N ) e^'^ , 

V K.N^ir.A = C w C K {Vhyp)^ (l + \t k ■ T N ) | ^000 ) (0000 1 J 

l^ooXCool . (7) 

where f are 1/2 of the Pauli isospin matrices for the 
1 = and 1=1 cases, and c v = \/E n f n (\/2TTa) 3 / 2 /\/3 
is a PCAC factor, and r is the relative coordinate. 

We calibrate our parameters in the two-body K • N 
channels, where the diagonalization of the finite differ- 
ence hamiltonian is straightforward. From baryon spec- 
troscopy we get (Vhyp) = 390MeV. Since ck is a 
PCAC factor that suppresses the norm of quasi-Golstone- 
bosons, consistent with the Adler zero in the chiral 
limit, we get the plausible parameter interval c|(Vft, yp ) € 
[200, 300] MeV. From the phase shifts for the first channel 
of eq. Q, t he repulsive and well understood S = 1 K»N 
channel|34. HH, [36| depicted in Fig. [3l we determine the 
size parameter a € [0.34, 0.39] fm. This is necessarily 
smaller, by a factor of 2, than the hadronic charge radius 
which is enhanced by the vector meson dominance. 



III. RESULTS AND CONCLUSION 

To study the binding and decay of the K~ • N two- 
body systems and of the K~ uNuN three-body systems, 
we diagonalize the hamiltonian in configuration space. 
Importantly, the decay widths to the channels n • E and 
7r • E • N are accounted with the substitution method, re- 
sulting in an effective two-body K~ • N separable poten- 
tial -|0ooo)(0ooo|G o (£)|0ooo)(0oool- Because tne energy 
E of our plausible resonances is above the pionic chan- 
nels, the Green function matrix elements are complex. 
This produces the decay width contribution — iY/2 to the 
eigenvalues of the Schrodinger two-body and three-body 
Hamiltonians. 

First we study the A channel with the 1 = two- 
body K~ • N system. Using the PCAC ratio c^/ck — 
\J E-k/Mk — 0.56, and solving the K~ • N Schrodinger 
equation in this parameter range, including the coupling 
to the 7r • E effective potential, we get a resonance in the 
A(1405) region, with binding energy M — ttlk — wjv £ 
[-23.8,-1.7] MeV, and width Y € [3.5,7.0] MeV. This 
is distant from the single Breit-Wigner model for the 
A(1405) resonance, with a width of Y ~ 40 MeV H3| 
and a binding energy of —30 MeV. We also consider 
an increase of the chiral coupling up to c^/cx = 1, 




FIG. 3: 1 = (grey) and 1=1 (black) s-wave K • N phase 
shifts as a function of the Kaon momentum in the laboratory 
frame. The dots are experimental l34l. [35l. [36| and the solid 
lines represent this model [29ll3CLl3ll.l32ll33| with the hadronic 
size parameter a = 0.36 fm and Ck 2 Vh yp = 436/ vo = 250 
MeV . A much larger a = 0.66 fm, even with a small 
Ck 2 Vh yp = 45 MeV potential strength, would produce a worse 
fit (dashed line). 



to simulate the possible effect of other decay channels. 
Maintaining the remaining parameter ranges, the reso- 
nance approaches the A(1405) region, with binding en- 
ergy M - m K - m N e [-42.2, -7.7] MeV, and width 
T € [11.2,23.0] MeV. This is reasonably close to the 
A(1405), possibly also close to its double pole models, 
including a slightly narrower Breit-Wigner pole closely 
below the K ~ • N threshold (here we don't address the 
second pole, since for simplicity we did not include the 
7r»S interaction, less relevant for the K~uNuN system). 

It is then instructive to find how much we can decrease 
the 1 = attraction, and still bind the K~ • N system. 
Excluding the coupled channels, the Vk-%n attraction is 
marginally sufficient to bind the 1 = system, up to a 1 
MeV binding energy The coupled channel effects further 
bind the system and, even if the overall 1 = attraction 
is reduced by a factor S [0.77,0.89], binding occurs. 

In the S channel with the 1 = 1 two-body K~ • N 
system, according to eq. ([7]), the Vk-,n attraction is 
decreased by a factor of |, when compared with the 1 = 
channel. Clearly, this small potential is far from being 
able to bind the K~ • N onto a E. 

We finally study the K~ • N • N three-body system, 
using the Nijmegen REID93 potential [H, together with 
our K~ • N interaction and our coupling to the -k • E and 
7r • A channels. Importantly, our two-body sub-systems 
only bind for isospin 1 = 0. For the N • N sub-system 
we can either have isospin Inn — and spin Snn = 1- 
or Inn = 1 and Snn = 0. The two different total I = \ 
systems are respectively the (— K~pn + K~np)/\/2 with 
spin 1, and the (— 2A > onn — K~ pn — K~ np) / \/Q with spin 
0. In the 5=1 case, each K~ • N sub-system is \ in 
a A-like 1 = state and § in a E-like 1 = 1 state. In 
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the 5 = case, each K • N sub-system is | in a A-like 
1 = state and \ in a E-like 7=1 state. For instance, 
in the extreme case where the two nucleons would be 
superposed, we can use eq. (J7]) to estimate the attraction 
that this dinucleon would impose on the antikaon. In the 
5 = 1 system, the K~ would feel an attraction identical 
to the one it feels in the A system. In the 5 = system, 
the K~ would feel a larger attraction by a factor of |. 

We now detail the study of the S = 1 K~»N»N, where 
the N»N sub-system binds in a deuteron. It is convenient 
to replace the coordinates of the three hadrons Ni, N2 
and K~a by relative and centre of mass coordinates 
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where 7 = mK/rriN- The coordinate R is eliminated in 
the centre of mass frame, and we rewrite the hamilto- 
nian in terms of scalar products of 172 , 1*123 and their 
momenta P12 and P123. The total kinetic energy, 

12 o 1 2+7 
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(9) 



is diagonal in P12 and P123. The Vn x n 2 potential only 
depends on 172, including the tensor interaction. The 
sum of the other potentials Vn 1 k 3 + Vn 2 k 3 is a function 
of ri2, 7*123 and of the square of u> = 772 • 7723. To span 
the possible quantum numbers of K~ • N • N system, we 
apply the total potential energy to a groundstate prod- 
uct of s-wave wavefunctions of 172 and of 1723. In what 
concerns the angular quantum numbers, the tensor part 
of the N • N interaction Vn 1 n 2 wm couple the s-wave in 
172 to a d-wave. The other two potentials Vn 1 k 3 + Vn 2 k 3 
couple the groundstate s-wave product to wave-functions 
with any even power of to. We address the binding of the 
three-body K~ • N • N system with two different appli- 
cations of the variational method. In the first method, 
we use a basis of Laguerre Polynomials in the radial vari- 
ables 772 and ?723, 
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adequate to study weakly bound systems, together with a 
basis of Legendre polynomials Pi(to) in the angular vari- 
able W. We use up to 10 excitations in each Laguerre 
basis and 5 excitations in the Legendre basis, sufficient 
to get for instance the deuteron binding. In the second 
method we use a finite difference method, with 40 excita- 
tions in each of the radial variables, and 3 in the angular 
variable. The binding and decay width of the K~ • N • N 
are determined from the lowest eigenvalue of the hamil- 
tonian, both in the polynomial basis and in the finite 
difference basis. 

It turns out that in the K~ • N • N, I = 1/2, 5=1 
channel there is no binding. The groundstate has bind- 
ing in the 772 coordinate, but no binding in the ri23 co- 
ordinate. In particular, the r\i part of the wavefunction 
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FIG. 4: Kaon energy, or real part of the lowest s-wave eigen- 
value of T\23 + T w + VN1K3 + Vn 2 k 3 , as a function of the 
nucleon-nucleon distance ri2 = |7\zVj — risr 2 \. This result is 
obtained in a large sphere with a ri23 radius of 7.5 fm. In the 
infinite sphere radius limit the energy essentially vanishes for 
ri2 > 2.0 fm. 



is localized and reproduces the deuteron wavefunction, 
while the ri23 part is extended over the whole size of the 
large box where we quantize the wave-function. In the 
limit where the size of the box is infinite, we get a bound 
deuteron pun and a free K~ . 

To study in more detail how far we are from binding, 
we consider the adiabatic limit of the 7=1/2, S = 
1 system, where the positions of the two nucleons are 
frozen, by fixing the 772 variable. Then the energy of the 
Kaon K~ can be computed, diagonalyzing in detail the 
terms in the hamiltonian depending on 7*123 an d w. The 
resulting smaller eigenvalue of Ti23 + T w + Vn 1 k 3 +Vn 2 k 3 
is shown in Fig. 31 using the parameters in our interval 
that provide the strongest possible binding. It occurs 
that the antikaon is bound to this frozen deuteron system 
only if 772 < 2.0 fm. 

Notice that this is much s maller than the deuteron 
double radius mean square \J (ri2 2 ) — 3.9 fm, which is 
quite large because at short distances the N • N subsys- 
tem suffers a strong repulsion and because the deuteron 
is weakly bound. Thus is would be very hard to contract 
the deuteron to a sufficiently small radius to bind the an- 
tikaon. Then, at these large distances, the antikaon es- 
sentially feels a double- well potential, as shown in Fig. [5] 
Importantly, such a widely separated double-well poten- 
tial only binds if any of the two wells is sufficiently deep 
to bind the antikaon. It occurs that this already happens 
if the K~ • N attraction in the 5=1 three-body system 
is arbitrarily increased by a factor <= [1.14,1.35]. Then 
we get a vanishing binding energy, and a decay width of 
T G [8.0,9.8] MeV. 

This can also be applied to the 5 = 0, K~ • N • N 
three-body system, where we have binding because the 
attraction in the K~ • N sub-systems is increased by a 
factor of 5/3 when compared with the 5=1, K~ • N • 
N three-body system. In particular we find a binding 
energy M - m K - 2m N G [-53.0,-14.2] MeV, and a 
decay width T G [13.6, 28.3] MeV to the tt • S • N and 
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FIG. 5: 3d perspective (colour online) of the antikaon wave- 
function, assuming two adiabatically frozen nucleons at the 
distance of ri2 = 2.5 fm, obtained with polar coordinates 
for the ri23, with finite differences, and with a radial lattice 
spacing of 0.1 fm. In this unbound case, the wavefunction is 
a deformation of the lowest positive energy harmonic. 

7r • A • N channels. The complex pole of this resonance 
is comparable to the one we get for the A(1405). 

To conclude we compute, starting at the quark level, 
the K~ • N interactions, constrained by chiral symmetry 



and by the crossing symmetry to the K + • N system. 
We find binding in the I — 0, K~ • N channel, possibly 
consistent with the double pole model of the A(1405), 
and no binding in the 1=1, K~ • N two-body E. We 
also study the I = 1/2, K~ • N • N three-body systems. 
In the 5=1, K~ • N • N system, we find no binding, 
although we almost get binding. We find a binding in 
the S = 0, K~ u N • N system, with a complex pole 
comparable to the one of the A(1405). Essentially, the 
main difference to the generally larger bindin g an d decay 
of very recent theoretical calculations [f| [Tq Il7l. [l9| 
can be linked to our shorter range and weaker potential 
Vk- n- 
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